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Abstract 

A new approach that uses non-structured mesh to replace the arbitrarily overlapped struc- 
tured regions of embedded grids is presented. The present methodology uses the Chimera com- 
posite overlapping mesh system so that the physical domain of the flowfield is subdivided into 
regions which can accommodate easily-generated grid for complex configuration. In addition, a 
Delaunay triangulation technique generates non-structured triangular mesh which wraps over the 
interconnecting region of embedded grids. It is designed that the present approach, termed 
DRAGON grid, has three important advantages: (1) eliminating some difficulties of the Chimera 
scheme, such as the orphan points and/or bad quality of interpolation stencils, (2) making grid 
communication in a fully conservative way, and (3) implementation into three dimensions is 
straightforward. 

A computer code based on a time accurate, finite volume, high resolution scheme for solving 
the compressible Navier-Stokes equations has been further developed to include both the Chimera 
overset grid and the non-structured mesh schemes. For steady state problems, the local time step- 
ping accelerates convergence based on a Courant-Friedrichs-Leury (CFL) number near the local 
stability limit Numerical tests on representative steady and unsteady supersonic inviscid flows 
with strong shock waves are demonstrated. 



I. Introduction 


During the last decade both structured and unstructured grid systems have been developed 
and applied for the computations in various CFD problems. Among various grid concepts, there 
are currently two mainstream approaches for computing flow fields in which geometry imposes 
complex boundary conditions - composite structured grid schemes and unstructured grid schemes. 
In general, unstructured grid methods are considered to be more versatile and easier to adapt to 
complex geometry while composite structured grid methods are considered to use more efficient 
numerical algorithms and require less computer memory. However, both pure-strain approaches 
have their strengths and weaknesses. 

The Chimera grid scheme 1 which uses overset grids to resolve complex geometry or flow fea- 
tures is generally classed into the composite structured grid catagory. Overset grids allow struc- 
tured grids to be used without excessive distortion or inefficient use of grid density. It has been 
used to compute inviscid and high Reynolds flow about complex configurations, 2 ' 3 and even dem- 
onstrated for unsteady three dimensional viscous flow problems'* in which one body moves with 
respect to another. Overset grids have also been used as a solution adaptation procedure. 5 * 6 It is 
found that the Chimera approach has more versatility than current unstructured schemes. 

Up to date, the Chimera method is an outgrowth of trying to generalize a powerful solution 
approach to more complex situations. However, it is also recognized that there are weakness 
which must be removed if Chimera is to remain competitive with unstructured grids. There are 
two main criticisms leveled against the current implementations of the Chimera method: (1) the 
complexity of the interconnectivity is perhaps as difficult as dealing with an unstructured grid, 
many times resulting in orphan points and/or bad quality of interpolation stencils left in the com- 
putational domain, and (2) nonconservative interpolations to update interface boundaries are often 
used in practical cases. The fact that interpolation is generally used to connect grids implies that 
conservation is not strictly enforced. For some practical applications, such as the turbine-blade 
cooling passage flow in which heat transfer is involved and the supersonic inlet design where 
shock waves travel in time, the interpolation errors become overall-significant if numerical fluxes 
are not fully conserved. 

It is remarked that while the space shuttle simulations have used simple interpolation proce- 
dures because of their robustness, Calspan/AEDC simulations have implemented an unpublished 
idea of Benek and use interpolants of delta quantities. Interpolating this quantity on interface 
boundaries ensures space-time conservation over the global field, but the penalty of this -approach 
has not yet been rigourously examined (or clearly addressed). Many conservative interface 
schemes for overset grids have been devised. 7 ' 8 More recently, Wang and Yang 9 transformed the 
conservation condition for overlapped grids to one of enforcing conservation for patched zones. 
However, a main drawback for those interface treatment methods is in determination of the area 
weights for irregular polygonal grid cells. Implementation into three dimensions is obviously dif- 
ficult and not straightforward. 

A fairly obvious way to ensure interface conservation would be to introduce an unstructured 
flow solver in the vicinity of the interface boundaries. Hybrid schemes which incorporate the best 
features of both structured and unstructured have already appeared 10 " 11 and will likely become 
more important in treating flow about complex geometries. In this paper a new approach is pre- 
sented which uses non-structured mesh to replace the overset region for Chimera embedded grids. 
This approach not only eliminates the above criticisms but also preseves the advantages of the 
Chime ra grid scheme. The extension of the present method to three dimensions is obviously 
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straightforward. 


II. Structured and Unstructured Grid Schemes 

Chimera Grid Scheme - Structured Grid System 

The Chimera scheme is a grid embedding technique which provides a conceptually simple 
method for domain decomposition. For instance, a major grid is generated about a main body ele- 
ment and minor grids are then overset on the major grid so as to resolve interesting features of the 
configuration. Usually, the minor grids are overset on top of the major grid without requiring the 
mesh boundaries to join in any special way. However, a common or overlap region is always 
required to provide the means of matching the solutions across boundary interfaces. To increase 
the flexibility in the selection of subdomains, the Chimera scheme also allows an implementation 
to remove regions of a mesh containing an embedded grid from that mesh. That is, an embedded 
mesh introduces a “hole” into the mesh in which it is embedded. Typically, a hole is defined 
through a creation boundary which consists of a surface or a group of surfaces. The purpose of a 
hole cre at ion boundary is to identify points that are within this boundary. A mesh point is consid- 
ered to be inside a hole creation boundary if it is inside all surfaces that define the boundary. Fig- 
ure 1 illustrates the connections between composite overlapping grids, with hole points being 
blanked by a prescribed creation boundary. Communication between overlaid grids is achieved by 
interpolation of boundary values from the mesh or meshes in which the boundaries are contained. 
For instance, the scheme employs nonconservative trilinear interpolation that as some simple 
experiments have shown is superior to Taylor series expansion. A practical application of the Chi- 
mera grid system for complete Space Shuttle/Solid Rocket Booster geometry is shown in Fig. 2. 
Details regarding the grid embedding technique can be found in Ref. 12. 



Figure 1. Interconnections of minor and major grids for Chimera overset grid scheme. 
(Ref. 12) 
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A flow solver must be modified to account for the use of multiple meshes and the hole in 
the grids. These hole points must be blanked or excluded from the flow field solution. The main 
change in the flow algorithm itself is the treatment of the hole boundaries. The hole information 
from the Chimera grid package is stored in an array, IBLANK, which is defined for each point on 

each grid as 


IBLANK = {*’ 


if a point is not blanked; 
if a point is blanked. 


( 1 ) 


In the flow solver, each element in the corrector step is multiplied by the appropriate IBLANK 
value before the solutions are updated: 


U n+ 1 = U n + IBLANK - (U n + l -U n ) 


( 2 ) 


The blanked solutions are updated in the interpolation routine. With this approach, no special rou- 
tine or logic tests are required to exclude the blanked points from the flow field solution. 



Figure 2. Overlapped grids for the Shuttle launch vehicle configuration. (Ref. 3) 


DRAGON Grid - Hybrid Grid System 

The application of the Chimera scheme involves two steps; (1) data communication between 
two embedded grids, and (2) execution of a flow solver that uses the communication information 
generated in step 1. The PEGSUS code is used to determine the interpolation coefficients between 
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composite grids. It is now modified to persue the DRAGON grid technique. The approach of the 
DRAGON grid is accomplished by using non-structured mesh in the vicinity of the interface 
boundaries to replace the arbitrarily overlapped structured regions. Descriptions regarding the 
procedures to adopt the DRAGON grid method are summerized as follows: 

(1) The entire domain of a minor grid is used to form the hole creation boundary. 

(2) Grid points inside of the hole creation boundary are blanked out in the major grid. 

(3) Overset grid meshes are no longer exist Instead, a new generated gap region is defined by 
the hole creation boundary as well as the hole boundary points. In result, interpolations of numer- 
ical fluxes or primitive quantities are no longer performed. (Fig. 3(a)) 

(4) Both fringe points (hole boundary points) and interpolation boundary points are no longer 
treated as blanked points. Instead, they are now represented as interior points, and their IBLANK 
values are saved to be 1 in the PEGSUS code. 

(5) The fringe points as well as the boundary points are now forming the new boundaries for 
the unstructured grid generator. Details of unstructured grid generation will be addressed later in 
this section. A file named DRAGON is used for output of PEGSUS code which contains the 
boundary mesh coordinates of the unstructured region. 

(6) Since interpolation process is not performed between structured grid blocks, output files 
providing interpolation information are deleted from the PEGSUS code. 

It is known that the triangular shaped cell is the simplest geometric shape to cover a computa- 
tional domain. A gap region created by arbitrary embedded grids implies an irregular shape of 
domain and would be very difficult to be represented by structured grids. The well-known 
Delaunay triangulation scheme 13 is applied to perform the unstructured grid generation. The steps 
which adopt the triangular cells in combination with Chimera grid scheme are summarized. 

(1) Boundary nodes provided by the PEGSUS code are reordered according to their geomet- 
ric coordinates. 


<•) 




Figure 3. Direct replacement of arbitrary grid-overlapping by non-structured grid, 
(a) Chimera grid method, (b) DRAGON grid method 
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(2) Delaunay triangulation method is then performed to connect the boundary nodes based on 
the Bowyer algorithm. 14 

(3) Triangles interior to bodies or exterior to the domain are deleted. 

(4) Interior nodes are inserted in the circumcenter of triangles which are deemed to be too 
large based on the boundary spacing and a clustering parameter (which governs the rate of change 
of triangle area). This step completes the specification of the initial grid. (Fig. 3(b)) 

At present, no further grid inserting or deleting methods have been performed. Mesh 
impr ovement techniques, such as grid smoothing, aspect ratio control, equal length sides, etc., can 
be used to achieve better grid quality. An example of a coarse grid system about the multi-body 
domain by using the Delaunay triangulation method is shown in Fig. 4. 



Figure 4. Delaunay triangular meshes for multi-body domain, (by Dr. James Loellbach) 


With an unstructured grid, individual cells can no longer refer to their neighboring cells sim- 
ply by using grid indices as in a structured grid. Instead, a connectivity matrix would be used to 
identify the relative position of the neighborhood of a cell. This connectivity matrix usually con- 
tains the cell based information as well as edge based information. With the present approach, 
connectivity matrices are also formed to define the discretized geometry between structured and 
unstructured grids. The connectivity matrices used in the present approach are summerized as fol- 


lowing: 

EEDGENODE(l:2, edges) 
ICELLEDGE( 1 :3, ceUs) 
IEDGECELL(1:2, edges) 
IEDGETYPE(edges) 
IEDGEFLUX(grids, edges) 


- 2 nodes of each edge 

- 3 edges of each triangular cell 

- 2 neighboring cells of each edge 

- edge type (type of boundary condition) 

- edge number that connects structured and unstru- 
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IFLUXINDXCgrids, edges) 
IFLUXINDY (grids, edges) 


ctured grids 

- /-index of structured cell that shares the edge fluxes 
with unstructured mesh 

-./-index of structured cell that shares the edge fluxes 
with unstructured mesh 


III. Governing Equations and Solver Algorithm 
Governing Equations 

Considering both structured and unstructured grid systems, the time-dependent compressible 
Euler equations can be expressed in an integral form in a volume v with enclosing surface s as fol 
low 


dv + §fids • F — 0 

V S 


(3) 


where the conservative variable U is 


V = [p,p2,p£] r 


(4) 


the inviscid flux F 


T 

F = [pu, puu+pl, pHii] 


(5) 


and the equation of state for ideal gas 


P = (Y-l)pe 


( 6 ) 


Here e represents the internal energy. The vector quantities, expressed in terms of Cartesian coor- 
dinates, are denoted with an overhead arrow , and the tensor with an overhead tilde or a dyadic 

notation such as uu. 

Flux Splitting 

Based on the finite volume method, the governing equation is semi-discretized by assuming 
that die cell-centered conserved variables are constant within a cell, and that the flux integral at 
cell surface is also approximated by an average value of the numerical flux and the surface length. 
The nume rical formulation uses a new class of flux splitting scheme. Accordingly, the name of the 
scheme is coined as Advection Upstream Splitting Method (AUSM). 15 Hie scheme first splits the 
full flux into convective and pressure fluxes, F\ /2 and F\ n respectively, 
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( 7 ) 


F\/2 ~ F l/2 + F l/2 ~ C U2^L/R +P t/2 

at the cell interface L<^<R. The interface pressure, P l/2 , simply comprises the positive and 
negative data from appropriate domains of dependence via characteristic speed decomposition. 
The passive scalar variables O s (p, p u, pH) is transported by a common convective velocity 
c that is constructed in a similar fashion as P l/2 . Then the upwind idea is used to select the 
state, i.e., “L” or U R”, of the variables to be convected. As such, the interface flux can be recast in 

the following form: 

-\\cuz\\n ts>+F uz (8) 

where A ()=()-() L . Here the first term on the right hand side is clearly not a simple 
average of the “L” and “/?” fluxes, but rather a weighted average via the convective velocity. 

Time Integration . 

An integration scheme assembles fluxes entering the cell interfaces and updates the conserva- 
tive variables at the cell center. The present method originates from the Taylor series expansion in 
time, as was done in the Lax-Wendroff scheme. A two-step scheme with second order time accu- 
racy is given as: 

predictor: 


U 


* 


U n + At- 


du n 

dt 


( 9 ) 


corrector: 


U = 




if 1+1 = \(U n + U**) 


( 10 ) 


The residual ^ vanishes as the solution approaches a steady state. This scheme is restricted by 
the numericarltabifity condition. Note that both predictor and corrector steps are identical, with 
no need of defining a midpoint for the corrector step. This reduces the complexity of evaluating 
die transport terms in the general curvilinear coordinates. 


Local Time Stepping _ ...... . ..... . . , 

For steady state cases, local time stepping accelerates convergence to steady state by advanc- 
ing the solution at each cell in time at a Courant-Friedricks-Leury (CFL) number near the local 
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stability limit. Thus, the expression for the local time step at each grid point (£, T), Q, 

At = CFLV/(\U]+a)S 


(ID 


where V is the cell volume, a is the local speed of sound, and S is the projected area in r\, and £ 
directions. 

Data Communication Through Grid Interface 

For multi-block grid system, the flow field is usually decomposed into separated regions 
according to different body geometry or flow features. However, the success of the overlaid grid 
or patched grid techniques relies on correct communication among grid blocks through boundary 
interfaces. Since the interface treatment methods are not necessarily satisfying any form of a con- 
servative constraint, the solutions on overlaid grids are often mismatched with each other. This 
may result in instability or oscillations, especially when a shock wave passes through boundaries 
of overlaid grids. 

In the current work, both the structured and unstructured flow solvers were based on a cell 
center scheme. The rectangles and triangles are the control volumes used in the finite volume for- 
mulation. Figure 5 illustrates the space discretizations of the volume cells for both structured and 
unstructured meshes. In the finite volume formulation (Eqn. 3), numerical fluxes will be evaluated 
at the cell interface based on the conditions of neighboring cells (L and R cells in Fig. 5). For 
unstructured grid, the edge flux ,/*, of the interface will be evaluated using the stuctured cell value 
as the right and the unstructured cell value as the left On the other hand, for structured grid, the 
inter fac e fluxes which have been evaluated in the unstructured process are applied directly in 
computation of the cell volume residuals, i & Urn j -ft- & * s noted that no flux interpolation has 
been performed at any interface over the entire multi-block domain. Obviously, the present 
DRAGON grid approach automatically preserves both the local and global conservation law, and 
further implementation into three dimensions would be straightforward. 



Structured/Unstructured 

Interface 


Figure 5. Control volumes for flux calculations. 
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IV. Tested Results and Discussion 


Computational results are presented in this section for two-dimensional, inviscid, steady and 
unsteady problems. Test cases include supersonic flow past a circular cylinder and the shock tube 
problem. Numerical solutions have been obtained based on the following grid systems: (1) single 
grid, (2) Chimera grid, and (3) DRAGON grid. Integration of the equations of motion uses a first- 
order-accurate version of the flux splitting scheme in both structured and unstructured regions. 

Case 1. Supersonic Flow over a Circular Cylinder 

The first problem solved was that of a circular cylinder in a supersonic free stream (M 00 = 2.5) 
with the associated bow shock. Three different mesh systems are employed. A single grid (Fig. 
6(a)) with 61 x 51 points is first used to provide an analytic solution. The flow field is also 
resolved using Chimera overset grid scheme while the major grid (41 x 101) covers entire 
domain and the minor grid (61 x 3 1 ) wraps around the cylinder body. Figure 6(b) shows the over- 
laid meshes with hole points being blanked in the major grid. Nonconservative bi-linear interpola- 
tion method is used at the grid interface which allows fluid properties to communicate between 
the major and minor grids. With the DRAGON grid approach, the blanked hole is extended 
beyond the outer boundary of the minor grid, as shown in Fig. 6(c), and a new gap region is 
formed between two embedded grids. Delaunay triangulation method is then performed which 
generates triangular meshes to fill in this region (Fig. 6(d)). The new generated unstructured mesh 
includes 390 nodes and 616 triangular cells. 



Figure 6. Computational grid systems for supersonic flow over circular cylinder. 

(M„ = 2.5) 

(a) single grid, (b) Chimera grid, (c) hole points blanked in DRAGON grid, 
(d) DRAGON grid 
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The values of the dependent variables at all the gnd points are set equal to their ffee-stream 
values initially. The equations of motion with different grid systems are integrated until the solu- 
tion converged to its steady- state value. Figure 7(a) shows the pressure contours obtained at con- 
vergence with the single grid mesh. Numerical solutions are also calculated based on both the 
Chimera grid and the DRAGON grid schemes, as displayed in Figs. 7(b) and 7(c) respectively, 
showing a comparable shock location to that predicted with the single grid system. The captured 
shock wave in Figs. 7(b) and 7(c) are seen to a slighter extent simply due to a coarser grid density. 
As noted in Fig. 7(b), both mismatched countour lines and pressure oscillations are observed in 
the vicinity of die grid interfaces. 

For the DRAGON grid approach, both structured and unstructured results can be displayed 
on the same plot using the Flow Analysis Software Toolkit 16 (FAST) visualization package on an 
IRIS workstation. As plotted in Fig. 7(c), the countour lines smoothly across the intersection 
boundaries; in fact, those lines seem to have a continuous slope across the interfaces. Note that 
small gaps interrupting the contour lines are observed along the interface boundaries as the cell 
center values are displayed in the structured grid region while unstructured grid values are ploted 
at their cell vertex locations. 



Figure 7. Pressure contours for supersonic flow over a circular cylinder. {P^^O.71 , 
Pmax-6-18, A P=0.28) 

(a) single grid, (b) Chimera grid, and (c) DRAGON grid. 
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Figure 8 compares the pressure distributions along the centerline of the flow field. It is shown 
that, for this test problem, both Chimera and DRAGON grid methods perform very accurate pre- 
diction for shock location and strength. 



Figure 8. Pressure distributions along the centerline of the flow field. — 2.5) 


Case 2. Shock Tube Problem 

It is of great concern that inaccurate shock speed and strength may be captured if the numeri- 
cal scheme is not fully conserved. An attempt is to test the time-dependent, plane moving shock 
problem on the embedded grid system. The shock wave is moving into a quiescent region in a 
constant-area channel with a designed shock speed M s = 4. Like the first test case, the computa- 
tional domain is composed using the single grid, the Chimera grid as well as the DRAGON grid 
schemes. As shown in Fig. 9(a), the single grid system containing a grid dimension of 801 x 21 is 
used to provide an analytic solution for shock location. Figure 9(b) illustrates the Chimera grid 
system which includes a major grid (upstream region) of size 41 x 21 , and a minor grid (down- 
stream region) of dimension 781x21. The grid lines between two overlaid grids are embedded in 
any arbitrary orientation. Upon the replacement of grid-overlapping by non-structured grid, the 
DRAGON grid for the moving shock problem is shown in Fig. 9(c). Note that all three grid sys- 
tems apply equivalent grid densities in the overall regions, except for that near the interface 

boundaries. . , , 

It is designed that the initial shock position is marked in the upstream region. (Fig. 9) The ini- 
tial conditions for the flow field are assumed to have the left and right values as 

P L = 4.57 
P L = 13.21 
u L = 3.125 
v L = 0 


P*= 1 
/>» = 0.71 

Ug = 0 

v R = 0 
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(a) 



t 


Initial Shock Position 


Figure 9. Computational grid systems for moving shock problem. (M s — 4) 
(a) single grid, (b) Chimera grid, and (c) DRAGON grid. 
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Figure 10. Pressure contours for moving shock problem. (500 iterations, At=0.1) 



Figure 11. Comparison of pressure distributions along the centeriine of the channel. 
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Figure 10 displays the pressure contours with various grid systems at a downstream location 
after the shock wave passing through the interface boundaries. It is hardly to identify any differ- 
ence on the contour lines between three grid systems. However, a critical comparison of the pres- 
sure distributions along the centerline of the channel, as plotted in Fig. 11, shows that the Chimera 
scheme predicts a faster moving shock in the tube. As a result, the present DRAGON grid 
approach accurately captures the shock speed which ensures the numerical solution is preserved 
in a fully conservative nature. 


V. Conclusion 


In this paper we presented a new approach, termed the DRAGON grid, that uses non-struc- 
tured mesh to replace the arbitrarily overlapped structured regions of embedded grids. It is 
designed to further enhance the flexibility for Chimera overlapping grid technique and to enforce 
conservative grid communication between embedded grids. Numerical results obtained for steady 
and unsteady problems indicated that the fluid properties have been accurately carried through the 
interface boundaries. Future plans include unstructured grid improving and 3D implementation of 
the DRAGON grid are now ongoing and will soon be reported. 
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